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We study in this paper the efficiency of different two-photon states of light to 
induce the simultaneous excitation of two atoms of different kinds when the sum of 
the energies of the two photons matches the sum of the energies of the two atomic 
transitions, while no photons are resonant with each individual transition. We find 
that entangled two-photon states produced by an atomic cascade are indeed capable 
of enhancing by a large factor the simultaneous excitation probability as compared to 
uncorrelated photons, as predicted some years ago by Muthukrishnan et al, but that 
several non-entangled, separable, correlated states, produced either by an atomic 
cascade or parametric down conversion, or even appropriate combinations of coher- 
ent states, have comparable efficiencies. We show that the key ingredient for the 
increase of simultaneous excitation probability is the presence of strong frequency 
anti-correlation and not time correlation nor time- frequency entanglement. 
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I. INTRODUCTION 

Quantum entanglement, and its inherent non local properties, are among the most fas- 
cinating and challenging features of the quantum world. In addition, entanglement plays a 
central role in quantum information [TH3]. Since its first description in the decade of 1930 
([6]), and in spite of the decisive contribution of J. Bell|7j and the subsequent experimental 
studies[8], entanglement stays even now as a rather mysterious and puzzling property of 
bipartite quantum objects. In particular distinguishing between effects related to genuine 
entanglement and those related to the quantum correlations measured on a single quantum 
observable is a difficult task|^9j, as can be seen for example by the great number of papers 
about quantum discord [inilll]. Some time ago, a paper was published [I2] which discusses 
how entangled states would be able to induce transitions in quantum systems that factorized 
states could not excite. The physical problem studied in that paper is therefore a good test 
bench to examine in detail in a simple situation the role of entanglement and of correlations 
not related to entanglement. This is the purpose of the present paper. 

The problem under consideration is the probability of two-photon two-atom (2P2A) ex- 
citation, in the situation where the two atoms are of different species and have different 
transition frequencies and the light to which the atoms are submitted is in general non- 
resonant for each one, but resonant for the system of two atoms. Two photon absorption by 
single atoms or molecules have been studied since 1931 [13] and remains a current subject of 
theoretical and experimental research [T3H18]. When the atoms have more than one interme- 
diate state many important features, including cross section cancellation and enhancement, 
are observed [15] . These features have recently been shown to be applicable in characterizing 
the quantum states of the absorbed two photon [TH] . 

It has also been known for a long time, that the two photons resonant excitation of 
two different atoms is indeed possible when the two atoms are interacting [19]. A nearly 
monochromatic light beam will have a resonant two photon absorption peak when tuned 
across the average frequency of the two atoms. Different experiments have since then con- 
firmed this theoretical prediction [201422] . In addition to a direct potential interaction be- 
tween the atoms, like the dipole-dipole, cooperative 2P2A has also been predicted for pairs 
of atoms inside an optical cavity [23] ■ In this case the physical interaction is mediated by 
the radiation background surrounding the atoms. Ref [12] addresses the case of two photon 
absorption in absence of interaction between the two atoms, with the excitation made using 
some particular entangled state of light. The conclusion of the authors of [12] is that in 
some situations entanglement can replace a real physical interaction, which is a far reaching 
statement and an important physical property related to entanglement. Surprisingly, this 
question did not attract much attention during several years. The same subject was also 
considered, but in the context of spin entanglement in Electron Paramagnetic Resonance 
by K. Salikhov|24j. More recently, [^ the related problem of interaction with pairs of 
broadband spectrum photons has been discussed. 



3 



In this paper, we will determine the probability of 2P2A excitation by different multi- 
modal states of light. From these results we draw conclusions on the respective role of 
entanglement and of correlations not related to entanglement in such a process. Section 
11 gives the general framework in which the problem is treated. Section 111 derives from 
a second order perturbation theory the relevant transition probability. Results for various 
different two photon states, introduced in section IV, are given in section V and VI. Finally, 
in section Vll we discuss different hypotheses for the physical origin of the enhancement of 
the 2P2A process. 

II. THE MODEL 

Let us first precise the model we are using and the notations. We consider two different 
two-level atoms labeled (1) and (2), having ground and excited states \gi) and \ei) {i = 1, 2), 
corresponding Bohr frequencies uji and spontaneous emission rates 7^ , interacting with a 
quantized field. We assume that that the transitions occur in times must shorter than the 
two atom lifetimes so that we can consider that the two excited states have infinite lifetimes 
(71 2 ~ ) and keep for ever their excitation. For the sake of simplicity we will assume that 
the light source is far from the atoms, so that the only non-empty modes are plane-wave 
modes having a single propagation direction Oz and a single polarization. In this situation 
one can use annihilation operators depending only on the frequency a{uje) — a^. At the 
rotating wave approximation the hamiltonian of the system is then given by : 

Hatom = f^ibibi + tvjj2h\h2, h = \gi) {ei\ 
t 

V ^hh\^ fi{uje)ai + hh\ ^ f2{(^e)a£ + h.c, (1) 

e e 

where fi{ooe) = —idi^/uj(J^2^£^SLe^^^'^/'^ — fiee^^^^^^'^, fie being a slowly varying function of 
the photon frequency, can be treated as a constant fi = fi{uji) ~ fiiojg). di is the electric 
dipole matrix element of atom i, S the transverse section of the beam which is focused on 
the atoms, Zi the position of atom i and L the length of the quantization box, the mode 
density in terms of frequencies Ue being 2nc/L. For simplicity, we will set ~ ^2 — and 
hence will not consider propagation effects. 

The evolution of whole system is described by a unitary operator U, p{t) = U{\gig2) {gig2\® 
Pq)U\ where po is the input light state, which can be either a pure two-photon state |^^) = 
Ylikq'^kq\^ki^q) or a mixed state in its spectral decomposition form po = Yli^iPtiY^ 
The probability of 2P2A excitation is given by : 

m = Y.P,{^,\{gig2\U\t)\e,e2){e,e2\U{t)\g,g2)\^,) (2) 
1^ 
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so that the transition probabihty is known when the field variable operator {gig2\U{t)\eie2) 
is determined and the incident light state po is known. The exact expression of the evolution 
operator is unfortunately not easy to obtain. To simplify our discussion, we will use lowest 
order perturbation theory, which is a good approximation for bi-photon states which do not 
carry much energy. 



III. SECOND ORDER PERTURBATION THEORY 



If we assume that the coupling between the light field and the two atoms is weak, the 
leading term in the evolution is f/^^^ = -h~'^e~'^^^°"'"'^"f^^^ dr dsV{T)V{s), where V 
is the coupling term in the interaction picture. One has 



(eie2|f/(^)|^7i^?2) = e-**(-^+-^+^//'') 5^a„a„^^„ (3) 



mn 



Amn = flioJm) f2iuJn) (4) 

SO that the leading term of the co-excitation probability ^ reads 

Pit) ^ E A*^^ 

mn TT{a]alamanPo) (5) 

jkmn 

In the case of a continuous frequency distribution of photons, one must replace the sum by 
an integral : 

{eie2\U^%ig2) = e"**^"^''"^"'''"^'^^ / [ dcOmdUr,aMa{coMmn (6) 

where a{u) in the annihilation operator of a monochromatic photon of frequency u. 

Note that the coefficient Amn is the product of two factors which represent the response 
of each atom to the field. When time t goes to infinity these two factors behave roughly like 
Dirac delta functions centered on individual atom resonances (we will precise this argument 
in section VI). The 2P2A excitation probability is indeed induced by the wings of the incident 
light spectrum which are resonant with the atoms. Consequently, if one photon is absorbed 
by one atom, there is no reason why the second photon should be absorbed by the second 
atom in a way correlated to the absorption of the first photon. In other words, the 2P2A 
excitation phenomenon has no reason a priori to have a resonant behavior when the 2P2A 
resonance condition wi + W2 = i^m + i^n is fulfilled. 

However, as the process is non-linear and involves two atoms, it can be enhanced by 
taking advantage of correlation effects between the atoms or between the photons : 

• A first possibility consists in introducing an interaction between the two atoms. Let 
Uu — uji + 002 be the maximal Bohr frequency of the two-atom system : if one photon 
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with frequency ojk is absorbed, then the two atoms will be more likely to absorb another 
photon with frequency Uu — w/c in a resonant two-photon process [l9l [27] . This was 
experimentally demonstrated in [2nH22] using nearly degenerate photon pairs. 

• A second possibility is to use correlated photons to interact with the two atoms. Let 
us consider a source that emits correlated photons : if a photon is absorbed by one 
atom, then the remaining atom will interact with its correlated photon with a higher 
probability, leading to enhanced 2P2A resonance (we will precise this argument and 
the kind of correlation needed in section VII). 

We will now precise these latter ideas by having a closer look at different possible light 
states likely to induce such a 2P2A transition. 

IV. ENTANGLED, CORRELATED-SEPARABLE AND FACTORIZED 

TWO-PHOTON STATES 

Before we go further, let us precise the different kinds of two-photon states that we will 
consider in the following. Starting from entangled pure quantum state having a density 
matrix po= |^)(^| of matrix elements pkWqq' = (Ifc, lq|po|lfc', Ig'); can construct others 
that have the same mean energy and the same single photon spectrum, and hence that would 
give the same transition probabilities for a single photon resonance. We choose two special 
cases that will allow a quantitative evaluation of the role of correlations : 

• The first one is defined as 

Pi = ^Trpfc Ollgpoinfe (8)ng = ^Pfcfcggllfc, lg)(lfc, Igl (7) 

k.q k,q 

where 11^ is the frequency projection operator 11^ = pi is the diagonal part of 

Po- It has lost any temporal field coherence and is time independent. It is actually a 
correlated- separable stoie[28], which results from the "disentanglement" of the previous 
one. It gives rise however to correlations between its two parties. 

• The second one is defined as 

P2 = TrpfcPoJIIfc Tl[UgPo]Ug = '^'^ Pkkq'q'\lk){h\®^ Pk'k'qq\lq){lq\ (8) 

k q k,q q' k' 

This is a fully factorized state, which does not give rise to any correlation whatsoever. 

These states will induce 2P2A excitation with respective probabilities Pi{t) and P2{t)- 
The two diagonal density matrices {pi, P2) describe two c.w. fields whilst the entangled 
pure state po describes a pulse, as a result, at a time t, the flow of energy having interacted 
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with atoms in each state is different. However, as shown in appendix B, this quantity is 
quite the same in each state when t = L/c. For comparison, we will take t = L/c through 
the whole paper. 

We are of course interested in cases where one observes an increase in the excitation 
probability, i.e. when P{L/c) ^ P2{L/c); if P{L/c) ^ Pi (L/c), then entanglement is indeed 
the key to efficient 2P2A transition, whereas if P{L/c) ~ Pi (L/c) ^ P2(L/c), correlations, 
of quantum or classical origin, are more important than entanglement in the present problem. 

V. 2 PHOTON 2 ATOM EXCITATION INDUCED BY DIFFERENT PHOTON 

QUANTUM FIELD STATES 

We will now examine the efficiency of various multimode light states for the simultaneous 
excitation of the two atoms. 



A. two quasi-monochromatic uncorrelated photons 

Let us begin by the simplest case : two uncorrelated photon wavepackets of mean frequen- 
cies LOa and Up, and respective spectral widths 7q, and 7/3 much bigger than the detecting 
atom spectral widths 71 and 72 , emitted by two uncorrelated atoms excited at the same 
time in the past and arriving at the detecting atoms position at t = 0, described therefore 
by the bi-photon state \ilj^^{t)) with 

It is the tensor product of two single-photon wave packets [SU] of duration 7"^ and 7^^. In 
the calculation of the probability amplitude in ([s]), we will replace the sum over modes by 
the double integral (L^/47r^c^) ff^d(xJkd(xJq, extend each integration domain to the whole 
real axis and use the residue theorem [30]. For the transition probability amplitude, one 
finds, when •jat ^ 1 and 7^3^ ^ 1 while keeping 7it ^ 1 and 72t ^ 1 : 

^11 _ L^fl{uJl)f2{uJ2)gaiuJl)gi3{uj2) .^Qx 

~ c2(u;i - + «7«)(w2 -(^13 + ilp) 
If we assume that the coefficients g do not vary with frequency the normalization of the 
two-photon state imposes that : 

2c^7^ 

gc^gp = — ^ — (11) 

so that the transition probability P^^ is 

pii ^ -Po7«7/3 /-,2\ 

where Pq = d\d2'jJiUJ2/h^£lc^S'^ = 367r^7i72C'*/cJia;|5'^. 
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1. double resonance 

Let us first assume that we are in the best possible situation, where the photons are 
separately resonant with the two atoms : Ua = oJi and = 102- The transition probability 
is then equal to 



which can be written in a more general way 



P'^n = ^ (14) 

where Sfr is the effective area of frequency distribution |cfcgp in the {uk,ujq) plane (see 
figures ([2])). This result turns out to be general and implies that all pure states having the 
same effective areas Sfr, entangled or not, will produce the same doubly resonant transition 



probability. Thus we regard (13) as a universal result under the double resonance condition, 
and its value will serve as a reference for all subsequent transition probabilities. 

In an actual experimental situation, one may take : 71,2 ~ IkHz, 7q,^^ ~ IMHz, S ~ 

4:7l'^C^/uJiUJ2, thus Pl)^ji ~ 97112 /^T^'^lalp ~ 10"'^. 



2. two-photon two-atom resonance 

Let us now turn to the 2P2A resonant case, where none of the two photons are resonant 
with the two atoms, but where the sum of their two energies almost matches the sum of the 



two atomic energies : cu^ + cu/s — Wi + ijJ2- The transition probability ( 12 ) has in this case no 
resonant variation as a function of the 2P2A detuning 6 = Ua + ojjj — coi — U2- When 6 = 
the transition probability is : 

pll _ Pplalp _ pll 7a7^ , . 

^2P2A ~ /\4: ~ /\4 ^ ^ 

where A is the smallest frequency mismatch between the emitting atoms frequencies and the 
detecting atoms frequencies, supposed to be much larger than the atomic widths. Without 
loss of generality, we have taken A = Iwq, — Wil = |a;2 — W/jj. 

We then conclude that the special case of 2P2A excitation probability by uncorrelated 
photons is also non zero for any couple of frequencies uJa,(jJi3, thus such a two-photon 
transition turns out not to be disallowed but simply induced by the wings of the two single 
photon frequency resonances. It is therefore very weak, as witnessed by the A~^ variation 
of probability. 

B. two photons produced by an atomic cascade 



Let us now envision the case considered in [12] of a two-photon light state produced by 
a three-level atom excited at a given time t = in the upper state that cascades down to 
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the ground state on two successive transitions of Bohr frequencies successively equal to Ua 
and w^. The corresponding spontaneous emission rates are 7a and 7/3. We assume that the 
emitted light is wholly directed in the Oz direction of atoms (1) and (2) (by means of a 
parabolic mirror for example). It is described by a bi-photon wave-packet with a coefficient 
c^g'^ equal, at a time t long compared to the lifetimes of the two transitions, to[T2l [30] : 

cT: = T S.Ma.iu;,) 

Here this entangled non stationary state is produced by a cascade, so that the photon of 
frequency Ug always arrives just after the photon of frequency Uk- In addition, the probability 
to have photons of frequency sum Uk + ujq close to Ua + (jJp is high. We have therefore an 
entangled state which is not only correlated in time but also anti-correlated in frequency. It 
is the time-energy analog of the position-momentum entangled state introduced by EPR, or 
of the field quadrature entangled state[31 EUl [HT] . 

Using the Residue Theorem, the transition probability amplitude reads exactly 



A' 



C2 UJp2-5 -i{ll3-la) 
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«7o 



(17) 



When 7]^ 2 ^ ^ ^ la\i ^^^^ decaying terms in Eq.(17) are negligible, leading to a 
compact expression 



A" 



s 



(18) 



where u 



is the frequency difference between frequency and frequency 



CO,. 



fi,^ = a, (3, 1, 2, k, q. 



1. double resonance 

Let us first consider here also the most favorable case, which is the double resonance 
(DR) situation. Keeping only the largest term, one obtains in this case for the probability 
amplitude when 7^^ ^ 1 and 7^t ^ 1 : 

L'^flMf2M9aM9pM 

c27„7^ 

Using the same assumption as in the previous calculation, one finds for the probability 

Pdr = ^ = Ph'n (20) 

It is time independent because we are considering times much longer than the two-photon 
pulse of duration 7"^ + 7^^. As it is equal to the probability obtained with uncorrelated 
photons, we conclude that entanglement does not help in the fully resonant case, but does 
not harm either. 
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2. two-photon two-atom resonance 
Let us now turn to the 2P2A resonance case. One obtains in this case for the probabihty 

Pn 



2P2A 



4c2 52 + 72 



(21) 



U2 - UJ(3 UJl - UJ/S 

This expression, already obtained in [SU] , shows that for this state the probabihty has in- 
deed a resonant character around the two-atom two-photon resonance 6 = 0. The transition 
probabihty P2P2A exact two-atom two-photon resonance is then : 

peas ^ Pq _ pll (r,cy\ 

''''^"7a7/3A2"^^^A2' ^^^^ 

One therefore finds that the transition probabihty is in the present case smaller than 
by a factor (7^/A)^ at exact 2P2A resonance, as expected because one is now less resonant 
than in the double resonance case. One finds more importantly that P2P2A ^■^ larg^i" than 
P2P2A' than in the two uncorrelated photon case, by a factor {A/jaY, which can be 
very large. This enhancement of the 2P2A transition probability is the main result of [12] 
: entanglement may indeed significantly enhance the two-photon two-atom process. To the 
best of our knowledge no experiment has been undertaken to show such a striking effect. 

It must be emphasized that the present considerations do not imply that the atom cascade 
entangled state is the only one likely to produce such a significant increase in the transition 
probability. This is the reason why we will now consider other light quantum states which 
may also be of interest in the present problem. 



C. Correlated and factorized states analogous to the atomic cascade 

Let us now consider the two states that have the same energy and the same spectrum 
that we have introduced in section (IV) , namely the correlated-separable state : 

= frVE r 2^ 2, 77 — ^^\2^ J h,i,){h,U (23) 
i;" K/3 + 7^) [(wg/3 + wfc«)' + 7^] 

and the factorized state : 

The first one corresponds to an atomic cascade for which the starting time is random, thereby 
averaging to zero all the off-diagonal time dependent terms in the density matrix, the second 
one characterizes a mixed state with two uncorrelated photons having the same spectrum 
than the initial cascade state. They give rise to the following transition probabilities : 
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At exact 2P2A resonance, we have Pi ~ P^^^7^c^t^/(A^L^) and P2 — PoRlaljila + 
7/3)c^t^/(A^L^). At any time t, one finds Pi ^ P25 since the spectral widths are much 
smaller than the 2P2A detuning. This fact shows that correlations play indeed an impor- 
tant role in the efficiency of the excitation. 

Note that Pi and P2 depend on time, as can be expected in a situation where the detecting 
atoms, which have an infinite lifetime, are submitted to a stationary quantum state, and 



therefore to c.w. light. In order to compare Pi and P2 to P|p|^ (equation (21)), which is 
induced by a pulse of light, we need to fix an interaction time t. It is shown in appendix 
B that the two atoms are submitted to the same energy flow at time t = L/c. One then 
obtains at this time and at exact resonance : 

p, ~ pii JjL ~ pca« . (07) 

We thus find the result that a correlated- separable state like pi can induce the 2P2A transition 
as efficiently as the entangled cascade state. This statement constitutes the main result of 
the present paper. 

Let us stress that pi, though not entangled, has indeed genuine quantum properties, 
being a mixture of single photon states which are highly non-classical. It displays strong 
correlations that we will study in more detail in section VI. 



D. two-photon state produced by parametric down conversion 

Let us now examine the two-photon state \^pdc) produced by non-degenerate parametric 
down conversion which has been under wide and in-depth investigation for many years. 
Because of its x*"^^ nonlinearity, a non-linear crystal submitted to a pulsed pump field of 
central frequency Ua+ujp and narrow bandwidth 0"^ emits a signal field (central frequency Ua) 
and an idler field (central frequency W/j). Let be the frequency width of the phase matching 
curve. For the sake of computational simplicity we will use a Gaussian approximation for 
both the laser lineshape and the phase matching curve. The crystal generates in such a case 
an entangled state which is described by a wavepacket with a coefficient [32] given by 

4f = Me e + le (28) 
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where M is normalized coefficient, satisfying 



27rc 



In the expression of (28), we have assumed that the pump laser pulse had a Gaussian 
temporal shape centered at time to ^ (^a^ + (^^^ to provide most of the photons a chance 
to interact with the two detecting atoms. The factor i in the second component originates 
from a relative phase (depends on the birefringence) which is set to be 7r/2 for the sake of 
simplicity in our case. 

Here we will also extend the double integral to the whole plane and find, when t is 
sufficient large [33]. the transition probability 



ppdc ^ ^p^ 



at, 



2a o ^2 
; e 



,2 I , ,2 
'^2a+"l,3 



+e 



(29) 



Let us also take into account the two mixed biphoton states {pi'^'^, p^'^) pertaining to the 



pure SPDC type II biphoton state (28) 



pf" = AT' 



2 J2 e ^ 



2 I 2 
?2 



|lfc, lg)(lfc, 1 



kq 



ndc ^ C 



El 



e ^ +e ^ 



\lk){U 




(30) 



(31) 



where ( = 1 + aj/ial + aj). The ffist one corresponds to a SPDC process in which all the 
off-diagonal time dependent terms in the density matrix are averaging to zero by random 
processes, while the second one characterizes a mixed state with two uncorrelated photons 
having the same spectrum than the initial SPDC state. When t is sufficient large, their 
corresponding transition probabilities read 



Pi 



vrPn 



2aj 



/3 



ppdc 
-'2 



TT- 



2 al 



+ e 




(32) 
(33) 



We will once again take t = L/c to he able to compare in a fair way the pulsed and c.w. 
excitations through the whole following discussions. 
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1. double resonance 

Let us first consider the DR situation with oja = uji and W/j = a;2- Keeping the largest 
term, one finds the probabihty 



P'dr = P^DR - ^Po ^ , . ; (34) 

Once again, we conclude that entanglement is not active in enhancing the transition prob- 
ability in the double resonance case. 

One also finds ~ when a a = 'Ja,o'i3 = 7^- In the following we will take 

this correspondences of spectral widths for comparisons. Henceforth, Pff^ or Pj^ji will be 
regarded as a reference in the discussions related to SPDC biphoton state. 



2. two-photon two-atom resonance 

Now we will turn to the 2P2A case. The transition probability P2P2A indeed a 
resonant character around 5 = 0. At the exact 2P2A resonance, it is equal to 

P'^pIa - i^A^«e-^^^/-^ (35) 

which is much smaller than for the atom cascade state because the factor A^/o"^ enters now 
as exponent in a Gaussian function and the detuning A is much greater than the spectral 
widths. 

For the factorized, uncorrelated state p^'^^ the transition probability in this case reads 

P'2tp2A - Ph'Ril + 2a^/a^)-^/^e-^^^^/-^ < e-^'^'/^-'^^'^PptAA (36) 

Thus, P2P2A is niuch greater than the probability given by the factorized state because of the 
scale factor e^^ /(°"q+°"/3)_ go we obtain in the Parametric Down Conversion configuration the 
same conclusion as the one drawn in [12] for the atomic cascade : the entangled state \^pdc^ 
is much more efficient for inducing a 2P2A resonance than the factorized, uncorrelated state. 
For the correlated-separable state p^^^ the transition probability reads 

^l,2P2A — ^DR^ — ^2P2A \y ' ) 

Thus, one has -Pf 2P2A ^ -^2 2P2yi- '^^^ same conclusion is found as in the cascade case : the 
correlated-separable state is as efficient as the entangled state to boost the 2P2A resonance. 
The fact that -Pf 2P2A much larger than Pf 2^2^ ^^"^ -^f 2'p2A — P2P2A once again shows 
that correlations, which are not necessarily related to entanglement, play indeed a crucial 
role in the efficiency of the excitation. 
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VI. ENHANCEMENT OF 2P2A RESONANCE FOR MORE GENERAL 

CLASSES OF LIGHT STATES 

We have so far studied interesting but specific states of light and showed an enhancement 
effect for some of them, entangled or correlated-separable. It would be interesting to consider 
now more general classes of light states. 



A. Light pulses starting at a given time 

Let us go back to the initial equations (|4]) and ([s]). They contain functions like [1 — 
exp{iuimt)]/ujim- As explained in the appendix A, even though this function does not act 
as a Dirac function when it is applied to integrations with any function, it indeed tends to 
2iiT6{ui — Urn) whcu t — > oo when applied to functions of Um that have a Fourier transform 
which is strictly zero for t < 0. Such will be the case here. 

The initial two-photon light state |\E') is the pure state : 

1^) = 5^Cfc,|lfc,g, (38) 

kq 

Let us assume that this state describes a " switched-on" light which is not vacuum only after 
time t = 0. One can then use the delta function approximation. The probability that the 
two atoms are found in the excited state at times long compared to the pulse duration is 
now 

P~ ^^|Ci2 + C2l|^ (39) 

Mathematically, if |ci2| ~ |c2i| this interference, which has been studied in the literature [TB] . 
may lead to strong variations according to the relative phase. According to the Cauchy- 
Schwatz inequality, one has 

0<P<2(^^(|ci2p + |c2in). (40) 
However, physically speaking, only one component between Ci2 and C21 dominates in the 



expression (39). This is because we have assumed that the quantities uJi,U2,uJa,Ui3 are 
sufficiently separated from each other but with a small 2P2A detuning 5 ~ 0, as a result, Ui 
should be closer to one of the central frequencies of the fields than to the rest one. Under 
this condition, one has 

p^^^i\c,,\'+\c,n (41) 

The correlated and factorized states pi,p2 analogous to the initial state |\E')(\E'|, give rise to 
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the following 2P2A transition probabilities : 

= Y^Mici2r+ic2ir) (42) 
p - 

P2 = -jt^ ^ {\cinCm2\^ + \c2nCmlf') (43) 
mn 

A discriminability index in the role of enhancing 2P2A transition probability is defined 
by the quotient Gp between P and Pi at t = L/c, 



n - P 



t=L/c 



\Cl2 + C21 


2 


|Cl2 


2 + 


C2I 


2 



(44) 



Thus, one finds < Gp < 2. The maximum value 2 is achieved when Ci2 = C21. 

One has Cp ~ 1 under the physical conditions we stated before. That is, the entangled 
and the correlated-separable state yield almost equal transition probabilities. This implies 
that the conclusion that we had drawn in the special previous cases is valid for a large 
class of two-photon states : correlated states are as efficient as entangled states in 2P2A 
co-excitation when they have delivered the same amount of energy to the two atoms. 

Another important discriminability index is the ratio between the two transition rates Pi 
and P2 : 

- ^ ''"''t.'!"'' ^ (45) 



The value of the enhancement factor G12 can be used as a witness for the correlation needed 
in such a problem. 

Note in addition that, while P is sensitive to possible destructive interference effects 
between C12 and C21, Pi is not. Therefore, the enhancement effect as indicated by G, and 
due to correlations not related to entanglement, turns out to be more "robust" than the one 
related to it. 



B. Coherent states 

So far we have only considered biphoton states of different shapes, which are all strongly 
non-classical objects, as they are produced by spontaneous emission or parametric fluores- 
cence which are specifically quantum processes with no classical equivalent. But one can 
also envision superpositions of two-mode coherent states of the form : 

l^coft) = ^Ckq\a{ujk)) ® (46) 

kq 

where \a{uk)) is the coherent state \a) in the mode of frequency cuk, a being the same 
complex number for all modes. 
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The calculation of the transition probability must be redone from the beginning. By 
using the approximation (a|0) ~ valid for |a| ^ 1, one finally finds : 

P,oh{t) = \a\^P{t) (47) 

where P{t) is the probability ([s]) obtained for two-photon states. Apart from the energy 
scaling factor the conclusions of the previous paragraphs hold in the present case, which 
looks much more classical than the previously studied ones, as such states can be produced 
by classical means. 



VII. WHAT KIND OF CORRELATION IS REQUIRED TO ENHANCE THE 

2P2A TRANSITION PROBABILITY ? 

We have found in the previous sections that the 2P2A transition probability depends 
crucially on the specific state of light used for the excitation, even when all the considered 
states have the same energy spectrum. The question we address now is the physical origin 
of an enhanced transition probability. We have seen that entangled and not entangled states 
may give comparable results, so a first answer to the question is obviously that entanglement 
is not at the origin of the effect, but rather some kind of correlation effect which is shared 
by entangled and not entangled states. 

Candidates likely to play a role in the present problem is time correlation and frequency 
correlation. We will now examine them successively 



A. Temporal correlation effect 



It is well characterized by the cross second order correlation function g^^ {t, r) 



T 



Tr[po^i-Hi)^iT^W]Tr[po4"H^)<^(^)] 



(48) 



Assuming that the amplitude of the single-photon electric field is a smooth function of Uk, 
one gets for the pure state |\E') = J2kq'^k,q\'^k, Iq), 

2 



kq 



iUf^T — iujqt 



(49) 



It is the two-time Fourier transform of the two-photon state. 



1. In the case of the cascade state (16|) 



n- 



9{T)9{t-T)e 



_ N -27cr-27^(t-r) 



(50) 
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(2) 

FIG. 1. Plots of the cross temporal correlation function [t^r). The left one is for the atom cascade 
biphoton state, with 7^, = 0.05MHz and 7/3 — 0.5MHz; the right plot is for the SPDC biphoton state, 
in which the pulse takes place around to = 30/is, and with = 0.05MHz, ap — 0.5MHz, oj^q, = 2MHz. 
Note the {t, r) asymmetry in the first figure, and fringes in the second one due to interferences from two 
temporal processes. In both plots, one finds significant temporal correlations along the diagonal line. In 
a real condition, the value of uj^a should be much greater, leading to a poorer graphic representation for 
interference patterns 

9{t) being the step function. We notice here a time asymmetry between t and r, 
expected in the case of a cascade in which the photon is always emitted before the 
ojp photon. 



2. For the SPDC state (|28|) : 



2 

d!' ^) = [1 + sin {t-T)\ exp 



(51) 



As can be seen on the figure ([I]), g\ (t, r) is in both cases significant only very close to the 
diagonal, which implies that both states exhibit strong temporal correlations, as expected. 
The width of the diagonal, which gives the characteristic time of this correlation, is equal 
to 7^^(o"^^) in both the cascade and SPDC cases. 



It is easy to see that for the correlated-separable states (23) and (30), there is no time 

(2) 

dependence for (7^ (t,r), and hence no temporal correlation, as expected from a c.w. time 
averaged state in which the photons arrive at any time. It is also the case for the coherent 



states (46). As these states give 2P2A transition probabihties comparable to the entangled 



state, we must conclude that the temporal correlation is not the physical origin of the 
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enhancement effect, nor the time ordering of the photons present in the cascade state. The 
physical reason is that, as we have neglected their spontaneous emission, the two detecting 
atoms have an infinite memory time, and hence they can be excited separately at any time. 



It is well characterized by the cross second order frequency correlation function g'^{uj,uj') 



equal in the pure state case to \c{ijJ,ijj')\'^ and to p{uj,uj') in the mixed state case. 

This quantity is plotted in figure ^ for the cascade and SPDC states, either entangled, 
correlated-separable or factorized. One observes that the frequency correlation functions 
take significant values only on the anti-diagonal for the left side plots, which implies that 
the corresponding states exhibit strong frequency anticorrelations. This is not the case for 
the right side plots. The width of the anti-diagonal, which gives the characteristic width of 
the frequency anticorrelation, is equal to 7o(cra) in both the cascade and SPDC entangled 
and correlated-separable cases. 

Let us note that the entangled cascade and SPDC states are the only ones in our list 
exhibiting simultaneously time correlations and frequency anti-correlations : one has in 
these states EPR-like correlations, revealed by a violation of the time-energy Heisenberg 
inequality [3T| [3l] when 7^/7^ (or aa/crp) <^ 1. 

The important point to notice is that such a frequency anti-correlation exists for all the 
states which exhibit 2P2A resonance enhancement, and is not present for the states which 
do not give rise to this effect. We are therefore led to the conclusion they the property needed 
to enhance the 2P2A excitation is precisely the presence of strong frequency anticorrelations 
in the quantum state. 

This conclusion, that we have demonstrated for the two specific examples that we have 
considered in the first sections of this paper, is far more general, as can be seen on the 
expression of the probability written for any switched on two-photon state. 

Equations (37), (38) and (39) show indeed that the probability of 2P2A excitation is 
proportional to the component of the density matrix of the two-photon state corresponding 
to the existence of one photon with frequency ui and one photon with frequency U2- This 
gives a simple interpretation of the problem : there is 2P2A excitation only when each photon 
of the two-photon state is resonant with the atomic transition of the atom it excites. This is 
expected, since we are considering that the atomic excited states have a very long lifetime, 
and therefore very narrow linewidths. Since the spectrum of each photon of the source 
has a much larger bandwidth, the probability of excitation is small. If the photons are not 
correlated in frequency, the probability of double excitation is proportional to the product of 



B. Frequency correlation effect 




Tr[po^i~H^')4'Hw)4"^^(^)^«^^(^')] 



(52) 



Tr[poii"^(u;')^i''H^0]Tr[po4'Hc^)4'"H^)] 
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FIG. 2. Plots of the cross frequency correlation function {oj]^,ijjq). a : entangled, correlated- 
separable and coherent cascade states ; b : factorized cascade state; c : entangled, correlated- 
separable and coherent SPDC states ; d : factorized SPDC state. In all plots 7a = cr^ = 
0.05MHz, 7/3 = ap = 0.5MHz, uj^ = l.bMRz^ujp = 3.5MHz. The color codes, in the unit of 
on the top left (right) are shared by a and c (b and d). The left side plots exhibit strong 
frequency anticorrelations along the line uok + ujq = Wq, +a;^, whilst in the right side plots, one finds 
no such a correlation. The type II SPDC biphoton source is non-degenerate and each photon has 
two distribution peaks, thus one sees two bright spots in the left side bottom plot and 4 bright 
spots in the factorized case in the right side bottom plot. In a real condition, the distances of the 
peaks in the bottom side plots are much greater, and the sizes of spots are much smaller. 



the probabilities that each photon has the corresponding transition frequency, and this yields 
a very small transition probability. But when the photons are anti-correlated in frequency 
such that the sum of their frequencies is equal to the sum of the transition frequencies of the 
atoms, when one photon is resonant with one atomic transition, the correlated photon will 
be automatically resonant with the other transition, and the probability of 2P2A transition 
will in general be much higher than in the non-correlated case. 
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We can say that the 2P2A transition occurs with a higher probabihty when the sum of the 
photon frequencies is found inside a small interval around the sum of the atomic transition 
frequencies, so that the enhancement is associated with the inverse of the variance of the 
|cfcgp distribution in the direction of the diagonal. 

VIII. CONCLUSION 

We have now elements of answer to the question raised in the title and in the introduction 
about the role of entanglement in the two-photon excitation process considered in this paper : 
We have shown that what is necessary for the enhancement of the transition probability is not 
precisely quantum entanglement nor temporal correlations, but frequency anticorrelation, 
which can be due to the presence of entanglement in the state, but also to correlations that 
are not related to entanglement. 

As any nonlinear process, like two-photon absorption in a single atom [35j, 2P2A transi- 
tion probability can be modified by changing the quantum state of light, and therefore the 
enhancement effect that we have studied in this paper is due to the partial optimization of 
the quantum state. 

We have not treated in this paper the important question of characterizing in a quanti- 
tative way the frequency correlation relevant to the present enhancement and relating it to 
its classical or quantum character though various quantum correlation witnesses such as the 
quantum discord. It will be addressed in a subsequent paper, together with the important 
question of the full optimization of the quantum state with respect to the 2P2A probability 
maximization, given a constant spectral energy distribution. 
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APPENDIX A : 

IS (exp(— zwt) — l)/{2i-KUj) 
A GOOD APPROXIMATION OF THE DELTA FUNCTION? 

Let us note Stioj) the function (exp(— ■ico't) — l)/{2mu). One can also write it as Stioj) = 
— smut/{27iijj) — cos cot)/ (27rc<;). Whereas the real part of St(co') is a sine function which 
tends indeed to a delta function when t — )• oo, the imaginary part, being not a peaked 
function whose area is constant, is not an approximation of the delta function. So in general 
St{uj) does not tend to the delta distribution when it acts on the general set of integrable 
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functions. However, it can be so on a smaller set of functions. This set includes for example 
all the odd functions in u, a subset which is not relevant for the present paper. We show in 
this appendix that St(a;) behaves also as a delta function when it acts on functions which 
have a Fourier transform which is strictly zero before t = 0. 

Let us consider a function F{t) that is zero for t < and admits a well-behaved Fourier 
transform f{uj). Then 

-| poo -| /"OO 

/(^) = ^J diF(i)e-* - — dtF(t)e-* (53) 

/oo 
dujf{uj)e-^^ (54) 
-oo 

where f{uj) is absolutely integrable, which excludes functions like l/{uj+i^) from the present 
discussion. Let us now calculate the integral 

da; ''''^^"''^^^ ~ V (^) - ^ fdr H diuf{uj) exp{-iujT) ^ i f dTe{T)F{r) (55) 



Then I ^ i drFij) = 27ri/(0) when t — ?> oo. This proves that St{uj) acts as a delta 
function for the set of functions that have a Fourier transform strictly null for t < 0. 



APPENDIX B : 

WHY DO WE TAKE t = L/c IN THE COMPARISON OF TRANSITION 

PROBABILITIES ? 

In order to compare the probabilities of transitions induced by pulsed and c.w. light in a 
fair way, we must be careful to take the same amount of energy flow J-'(t) on the detecting 
atoms in both cases. This quantity is nothing else than the integral over time and transverse 
section S of the Poynting vector. It is equal to, at a given time t and for a state p : 

F{t) = 2eQcS [ Tr[pE+''{r)E+{r)]dr ~ Hujy [ Tr[pb^ {r)b{r)]dr (56) 
Jo L Jq 

where 6(t) = exp(— ia;,„T) and cu is the mean frequency of the state under consider- 

ation. 

For any diagonal density matrix(DDM), since Tr[p/)£ij\^6'''(r)6(r)] = 2 is time- independent, 
one finds a linear relationship between the energy flow and time t 

J^DDuit) = 2/10;^ (57) 

For any entangled pure state |^) = ^kq^kq\^k, Ig) '■ 

Tr[|^)(*|6t(r)6(r)] = J] | J^c^.e-'^^ + | 5^c,,e— ^ (58) 
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The energy flow at time t is 



drTr[|^)(^|6t(r)6(r)] 



(59) 



when most photons arrive at the detecting atoms after t = 0. One assumes that at suflicient 
large time t (much greater than the temporal coherence length of the fleld), the photons in 
state |\E') have fully interacted with the detecting atoms, therefore, one extends t to +oo 
without introducing notable error. By using the Parseval identity, one has 



J'^(t) ^ fujj- 

= hu- 
= hu) 



cirTr[|^)(^|6t(r)6(r)] 



— oo 
oo 



E 

, kq 



dr 



^kq 



E 

kq 



^kq 



2hhj 



(60) 



as expected. By comparison with Eq.(57), one flnds that at time t = L/c, the energies 
supplied by the c.w. fleld and by the pulse are equal. Under this situation, one can make 
reasonable comparisons between the corresponding transition probabilities. 
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